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We present a self-consistent theory, as well as an illustrative application to a realistic system, of 

phase control of photoabsorption in an optically dense medium. We demonstrate that, when propa- 

■j-»r ' gation effects are taken into consideration, the impact on phase control is significant. Independently 

r--s I of the value of the initial phase difference between the two fields, over a short scaled distance of 

f^ ■ propagation, the medium tends to settle the relative phase so that it cancels the atomic excitation. 

^SJ ' In addition, we find some rather unusual behavior for an optically thin layer. 

PACS number(s): 32.80.Qk, 32.80.Rm, 42.50.Gy, 42.50.Hz 

jKf-\ _ After the initial ideas ||l| and experimental demonstration |g,|[ of the feasibility of the control of photoabsorption 

and its products through the control of the relative phase of two fields, much work in atoms B and molecules 0] 

^-H ■ has explored a variety of processes. Many interesting issues &0] have been raised and clarified, establishing thus 

^ ■ the idea as a useful tool. With the exception of one paper H, however, theory and experiment have dealt only with 
single-atom (molecule) situations. But, if the ideas are to be contemplated for applications, the issue of propagation 
f^ . is crucial. Addressing this issue is the purpose of this Letter. 
[~^ ■ We have chosen the fundamental scheme MM which has served as a benchmark for many of the initial and continuing 

^^ r work. We consider the excitation of a bound transition through the combined effect of a single- and a three-photon 
^^ ' transition via two fields whose relative phase is controlled externally. We formulate and examine the propagation of 
<^ r a bichromatic electromagnetic field E through an optically dense medium consisting of Xe atoms. This electric field 
^S^ is a function of time t and space coordinate z and is composed of the fundamental and its third harmonic fields that 

SpH. have the same (linear) polarization and frequencies tuf and Uh — Sojf, respectively. It is expressed as 



o 



§ ■ E{z,t) = l[£;^e*('=/^-'^/*) +£;^e*('=''^-'^"*) -He. c], (1) 

O^" where Ej = £je *'^J , j — f,h, with £j and (pj the slowly varying in time and space real amplitude and phase of the 

^ I corresponding field, and kj = lujUjc^^, with rij the refraction index of the host medium at frequency ujj. Although 

T '~j ' in our present treatment the host medium is vacuum, and thus rif = n^ = 1 and kh = 3k f , for the sake of generality, 

rS I e.g., presence of a buffer gas, we shall keep in the formalism the refraction index. The frequencies uJh,f are chosen 

j^ ■ so that one harmonic photon and three fundamental photons are at near resonance with the transition from the 

ground state ( |1)) to the 6s state ( |2)) of Xe. A two-photon transition due to the strong fundamental or one-photon 

transition due to the harmonic fields lead to the ionization continuum (states |c)) of the atom. As we intend to 

explore intensities of the fields for which the one- and three-photon transition amplitudes between states |1) and |2) 

are of comparable magnitude so as to maximize the modulation depth, the transition |2) — > |c) would be dominated 

by the two-photon process and the one-photon ionization due to harmonic photon can be neglected. Experimental 

contexts for the situation we are considering have been detailed in the literature. 

Beginning with the second order wave equation for the field E{z,t), in the slowly varying (during an optical 
cycle) amplitude approximation one neglects all second derivatives and, after projecting onto the corresponding mode 
function exp[i{ujjt — kjz)], j ~ /, h, one arrives at 
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where Pj — P'e ^"^^ (P' being complex) is the slowly varying in time and space field-induced medium polarization at 
frequency uij. Consistently with Eq. (p, it can be expressed as 

P{z,t) = l[p^e'^kfz-^ft) ^ p^^^ik^z-u.^t) _^^^ ^j^ (3) 



The most general approach to the calculation of the response of the medium is through the atomic density matrix p 
which obeys the equation dtp = —ih~ [H^tom+D, p], with ifatom the free atomic Hamiltonian and D = —pE the atom- 
field interaction in the dipole approximation where p is the electric dipole operator. Introducing the rotating wave 
approximation and adiabatically eliminating the continuum and all virtual (nonresonant) bound states (connecting 
by the lowest-order paths the states |1) and |2)), the slowly varying density matrix elements of the two remaining 
states (Til — Pii, 1722 — P22, and cr2i — P21 exp[J3(cj/t + (j)f — kfz)], are found to obey the following set of equations: 
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(X {P2c If)'^ is the 2-photon ionization rate of |2) being proportional 



where 7 is the radiative decay rate of level |2), 71 

to the square of the intensity // ex £? of the fundamental (p2c is the effective 2-photon matrix element for the 

fundamental field on the transition |2) — + |c)), p\2 is the effective 3-photon matrix element for the fundamental field 
on the transition |1) — > |2), si and S2 are the lowest-order Stark shift coefhcients (polarizabilitics) of levels |1) and 
1 2), respectively, and pi2 = {l\p\2) is the matrix element of the electric dipole operator p. Finally A is the detuning 
of both fields from the |1) — > |2) transition resonance and 6 = {(f>h — 30/) — {kh — 3A;/)z their relative phase. 

Consider now the polarization P{z, t) = NTi[pp] of a medium of atomic density A'^. In expanding the trace of this 
equation, we again follow the same procedure as in obtaining Eqs. (0), i.e. we use the adiabatic approximation to 
expresses all density matrix elements that do not refer to the states |1) and |2) in terms of the three main elements 
cii, (722, and (J21- Equating the result with Eq. (|^), identifying and grouping together terms oscillating with the same 
frequencies, we obtain 



P; = 2N[£f{siaii H- 52^22) + ?,p''i^£}(J2i 
Pi = 2Npi202ie''. 
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Those equations, together with the Maxwell's Eq. (0) and the atomic density matrix Eqs. (0), provide a complete 
description of our system in terms of a closed set of equations. 

To present the numerical results for Xe, we use the parameters calculated previously [p[ via MQDT and appropri- 
ately converted to conform to the present definitions. For illustration purposes, it is desirable to have a maximally 
pronounced interference of the fundamental and harmonic fields. The respective Rabi frequencies are given by the 
first and second terms in the parentheses of Eq. (Ha). To obtain, for example complete cancellation at 6* = tt when 
these two terms are purely real and have opposite signs, it is obvious that the peak values and the temporal widths of 
both Rabi frequencies should be equal so as to overlap completely. Let the strong fundamental field have a Gaussian 
temporal profile with a peak amplitude E'^'^^ = £f[t = imax) and width tj. Then the peak amplitude and width of 
the weak harmonic field should satisfy the relations 
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In Fig. |l| we plot the ion yield Q — [1 — aii{t) — o'22(i)]t^oo at z = as a function of the relative phase 6 for 
three different intensities // of the fundamental. In all cases, the detuning A is taken such that it compensates the 
relative Stark shift of levels |1) and |2) at the maximum t^ax of the pulse, the harmonic pulse duration r/i = 1 
ns, and the conditions (|6|) are satisfied. In this figure, for all intensities and relative phase 9 = tt, the ionization 
vanishes completely since the two transition amplitudes interfere destructively and the second term on the rhs of Eq. 
(pa|), responsible for the stimulated transition from |I) to |2) is equal to zero throughout the duration of the pulses. 
Consequently, the medium practically does not interact with the fields and the atoms are "trapped" in their ground 
state |1). The more surprising result, however, is that in the case I'?'^^ = 8 x 10^° W/cm^, maximal ionization is 



found not for 9 = 0, 27r, as one would expect and is the case for the other intensities. This is a manifestation of 
the quantum-mechanical interference resulting from the fact that for this set of parameters, in Eqs. (Q) the terms 
responsible for the stimulated transition reach the maxima at ~ tt ± 0.287r where the ionization peaks are located. 
The numerical simulations also show that, while keeping the conditions m) satisfied, with decreasing pulse duration 
Tf, the ion yield reduces and its peaks a.t 9 ^ 0,27r, ... gradually disappear, which is analogous to the decreasing of 
intensity since the total energy of the pulse lessens. Increasing the intensity, however, results in a narrower dip in the 
ionization profile and a shift of its peaks towards the values of 9 that are closer to tt. 

Let us turn now to the propagation effects. The results presented below are obtained for a density of atoms N = 10^'^ 
cm^''. This, however, does not imply any limitation on the generality of the discussion since, as one can easily verify, 
the parameter ziVS, where S is the laser beam cross-section, is a propagation constant, and thus it is always possible 
to rescale the problem to any desired density and propagation length z. Conditions (m are assumed at the entrance 
to the medium. As we have noted above, in the case of initial phase difference 9{0, t) = tt, the atoms stay in the 
ground state and the medium appears to be "transparent" to both fields; neither the fundamental, nor the harmonic 
experience any remarkable distortion of their shapes or total energy Sj{z) oc J fii|£j(z,t)p, j = f,h, over distances of 
propagation z as large as ~ 50 cm. The accumulated over this distance change of the relative phase is only ~ IO^'^tt 
rad, which is due to the field independent phase shift of the fundamental, given by the term in parentheses of Eq. 

Consider next the case 0(0, t) = 0, i.e., at the entrance to the cell the two fields interfere constructively. The results 
corresponding to the parameters of Fig. |l| with J?'^'^ = 8 x 10^° W/cm^ are collected in Figs, y and |^. One can see in 
Fig. ^ that, in the course of propagation, the relative phase 9 (taken at the dynamic pulse maximum imax + -2/c) grows 
rapidly and over a distance of the order of 1 cm reaches the value tt, at which the initial constructive interference of 
the two fields turns to destructive. At the same time, the total energy of the harmonic pulse, after a small reduction 
over a short interval of z, begins to increase as a result of the energy transfer from the strong fundamental field, in 
the parametric conversion process. This small reduction of the harmonic takes place only at the beginning of the 
propagation, when the relative phase is still close to and the two fields interfere constructively, in the process of 
excitation of atoms from the ground state |1) to the state |2), while the generated part of the harmonic field is out of 
phase with the fundamental approximately by tt and continues to build up with a slight oscillation around the value 
TT of the phase. It is important to mention that throughout the propagation, the amplitude and the phase of the 
fundamental field do not change significantly. This is because the number of photons contained in that pulse exceeds 
by many (> 6) orders of magnitude the number of atoms the pulse interacts with over the distance of z < 20 cm. 
Comparing the three graphs of Fig. 0, one can see that with increasing of 9 and Sh, the ionization probability first 
also grows, which is consistent with the previous discussion related to that intensity of the fundamental filed. But as 9 
approaches tt, the ion yield drops almost exponentially until Q ~ 10%. This residual ionization that is present even at 
6* ~ TT (and tends to rather slowly) is caused by the fact that, because of the significant increase of the total energy 
of harmonic field, conditions (g) are not completely satisfied and the upper atomic level |2) acquires population due 
to that fraction of the generated field which exceeds the initial. Since the temporal widths of the pulses are less than 
the (radiative) relaxation time of the atomic coherence tT2i (7~^ — 2 ns), a significant fraction of the harmonic pulse 
amplitude is generated behind the fundamental (Fig. y). That part of the amplitude is then attenuated due to the 
atomic relaxation. Thus the total energy of the harmonic, after passing a maximum at z ~ 5 — 7 cm, decays then 
slowly back. Under these conditions, the leading part of the harmonic pulse that falls under the temporal shape of 
the fundamental is by ~ tt out of phase with the latter and therefore the ionization vanishes, while the generated 
tail is continuously scattered by the atoms in the process of radiative decay. The oscillations of the relative phase 
around tt are also slowly damped and the propagation reaches a "dynamic equilibrium" 

We note finally that a similar behavior of the system is obtained for a range of intensities we have explored. The 
main difference is that for weaker fields (/'?'*'' = 3 x 10^" and 1 x 10^" W/cm'^) the ion yield does not exhibit a 
maximum other than at z — and drops to zero much faster as z increases, which is consistent with the discussion 
above. 

We have examined the problem of propagation in the simplest context of phase control, namely the excitation of 
a bound state, which has been used as a prototype in much of the initial work |3j3]. As discussed here, the problem 
bears resemblance to earlier works ItMioI on cancellation in third harmonic generation experiments, and so does the 
whole issue of phase control. The relevance and possible impact of propagation has been recognized by Chen and 
Elliott B who presented data and an interpretation in terms of rate equations ||lQ]. Their study showed evidence 
of non-linear coupling, such as those discussed above, and called for "more rigorous techniques" in the approach to 
this basic problem. In the limit of validity of rate equations, our results do indeed recapture the equations employed 
in their analysis. It will be interesting to explore this issue under more general conditions, such as the excitation 
of states embedded in continua, on which we expect to report elsewhere. The basic features of our analysis should. 



however, remain valid. 

In summary, we have shown that the propagation of a bichromatic field with a preselected initial relative phase, has 
a profound effect. Over a rather short scaled distance and independent of its initial value, the relative phase settles 
to a value that makes the medium transparent to the radiation, precluding thus further excitation and consequently 
control. The scaled distance ziVS does of course involve the density of the species and the cross-section of the laser 
beam, which suggests some flexibility on the choice of these parameters. In any case, however, the actual length of 
the interaction region over which control can be active will be defined and limited by the combination of the above 
parameters, as well as by the geometry of the focused or unfocused laser beam. Briefly, for not very low atomic 
densities {N > 10^^ cm^'^), the harmonic field settles to the steady-state value within a thin layer where a focused 
beam is well approximated by a plane wave. In the presence of large ac Stark shifts, however, a detailed analysis 
including specific experimental parameters is mandatory. 
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FIG. 1. Ion yield Q = (1 -an - a22)t^oo ^ 1 ~(7ii{t 
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FIG. 2. Relative phase 9{z,t = fmax + ^/c) (a), normalized energy Shiz)/ShiO) of harmonic field (b), and ion yield Q{z) (c) 
versus propagation length z for the case I^^^ — 8 x 10^" W/cm^. 
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FIG. 3. Temporal profile of the amplitude £h of harmonic field at different z. All parameters are as in Fig. H. 



